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ABSTRACT 

In this paper, we develop a theory of redshift distortion of the galaxy power 
spectrum in the discrete wavelet transform (DWT) representation. Because the 
DWT power spectrum is dependent of both the scale and shape (configuration) 
of the decomposition modes, it is sensitive to distortion of shape of the field. On 
the other hand, the redshift distortion causes a shape distortion of distributions 
in real space with respect to redshift space. Therefore, the shape-dependent 
DWT power spectrum is useful to detect the effect of redshift distortion. We 
first established the mapping between the DWT power spectra in redshift 
and real space. The mapping depends on the redshift distortion effects of (1) 
bulk velocity, (2) selection function and (3) pairwise peculiar velocity. We 
then proposed /3-estimators using the DWT off-diagonal power spectra. These 
/9-estimators are model-free even when the non-linear redshift distortion effect 
is not negligible. Moreover, these estimators do not rely on the assumption 
of whether the pairwise velocity dispersion being scale-dependent. The tests 
with N-body simulation samples show that the proposed /^-estimators can 
yield reliable measurements of (3 with about 20% uncertainty for all popular 
dark matter models. We also develop an algorithm for reconstruction of the 
power spectrum in real space from the redshift distorted power spectrum. The 
numerical test also shows that the real power spectrum can be well recovered 
from the redshift distorted power spectrum. 

Subject headings: cosmology: theory - large-scale structure of universe 
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1. Introduction 

In three previous papers we have developed a method of measuring galaxy power 
spectrum with discrete wavelet transform (DWT) decomposition, which is an alternative 
of the Fourier power spectrum detection (Fang & Feng, 2000 (paper I), Yang et al 2001a 
(paper II), 2001b (paper III)). The DWT power spectrum estimator is information-loseless, 
and optimized in the sense that the spatial resolution is adaptive automatically to the scale 
to be studied. A test with observed sample of the LCRS galaxies showed that the DWT 
estimator can give a robust measurement of the power spectrum. 

In this paper, we continue our effort in this direction. The topic this time is to develop 
theory and algorithm of redshift distortion in the DWT representation. 

In terms of power spectrum measurement, the central problem of redshift distortion 
is to find the mapping between the power spectra in redshift and physical spaces. In our 
paper II, we have already studied this mapping. However, that mapping is directly obtained 
by a wavelet transform of the mapping of the power spectrum in the Fourier representation 
(Peacock & Dodds 1994). Although a Fourier mode can be transformed into DWT modes 
and vice versa, we should be careful in doing the transform of a second order statistical 
quantity from the Fourier representation into the DWT one. As we have showed in paper I, 
the covariance of density contrast in the DWT representation (ej^ej/j/) is not equivalent to 
the Fourier counterpart (^k^k), because the Fourier mode is subjected to the central limit 
theorem, while mode localized in both scale and physical space may not be so. This is 
because the DWT mode is characterized by not only the scale, but also the configuration 
(or shape) of the mode. Simply speaking, the variance (ej,iej',i') contains information of 
both the scale and phase (shape), while (^k(5k) does not contain information of phases (Fang 
& Thews 1998). 

We have found in paper II that, the Fourier mapping of redshift distortion can be 
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employed for diagonal DWT power spectrum, but not off-diagonal DWT power spectrum. 
This is due to the shape (phase) -dependence of DWT modes. A 3-D Fourier mode with 
wave-vector k = [ki, k2, k^) can be transformed into mode k' = [k'l, /cg, k'^) by a coordinate 
rotation as long as kf + k2 + k^ — /cf -|- k'2 + kf. Therefore, for an isotropic random field, 
the Fourier modes with the same k — |k| are statistically equivalent. However, the DWT 
modes do not share the same property. The length scale of a 3-D wavelet mode with scale 
index j = (ji, J2, J3) is L/[2^-'i -|- 2'^^^ + 2^-'3]^/^, where L is the length scale of the sample. 
Generally, one cannot transform a mode (ji, J2, J3) to (ii, J2) J3) by a rotation, even if they 
have the same scale. 

As a consequence, the redshift distortion in the DWT presentation will be shape- 
dependent. For instance, in the Fourier representation, the linear redshift distortion factor 
(1 -|- jS/i'^y (Kaiser 1987), where /3 is the so-called redshift distortion parameter, depends 
only on fjL — ks/k, i.e. the cosine of the angle between the wavevector k and the hne-of-sight. 
While in the DWT representation, the linear redshift distortion factor will depend on not 
only the scale, but also the shape of DWT modes. 

An important application of redshift distortion is to determine the redshift distortion 
parameter f3 (e.g. Hamilton 1998), which contains valuable information of the cosmological 
mass density parameter and the bias of galaxies. The shape-dependence of redshift 
distortion is very useful for the parameter determination. We will develop an algorithm of 
P estimation with diagonal and off-diagonal DWT power spectra. 

Moreover, the DWT representation provides an easy way to study the effect of 
selection function on the redshift distortion. In the first three papers, we assumed that 
the selection function n(x) is known. Actually, the "known" is not necessary. According 
to the definition, selection function is an observed galaxy distribution if galaxy clustering 
is absent. Therefore, for a consistent algorithm of power spectrum, the selection function 
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should be and can be obtained by the galaxy distribution itself. We will show that in the 
DWT representation the selection function can be determined by the distribution of galaxies 
without other assumption. This algorithm is convenient to estimate the contribution of 
selection function to redshift distortion. 

The paper will be organized as follows. §2 is a summary of the algorithm of the DWT 
power spectrum. In §3, we develop a theory of the redshift distortion with a multiresolution 
analysis. The emphases are the real-redshift mapping of diagonal and off-diagonal power 
spectrum. We proposed /3-estimators, which are tested by N-body simulation samples of 
popular dark matter models (§4 and 5). Finally, the conclusions and discussions will be 
presented in §6. The mathematical stuffs with the relevant DWT quantities are given in 
Appendix. 



2. Algorithm of the DWT power spectrum 

2.1. The DWT power spectrum 

We summarize the algorithm of the DWT power spectrum. The details can be found 
in Fang & Feng (2000) and Yang et al (2001). 

If the position measurement is perfectly precise, the number density distribution of 
galaxies can be written as 

n^(x) = J^'^iSoix - X,) = n(x)[l + 5(x)] (1) 

i=l 

where Ng is the total number of galaxies, Xi is position of ith galaxy, Wj is its weight, and 
Sd the 3-D Dirac 5 function. n{x) is selection function given by the mean number density 
of galaxies when galaxy clustering is absent, and 5(x) is the density contrast fluctuation in 
the underlying matter distribution. 
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Without loss of generality, we enclose the sample in a cubic box with side lengths 
(Li,L2,I/3). In the DWT representation, n^(x) is decomposed into 

7 oo oo oo 2^1-12^2-124-1 

n^(x)=n(j)(x) + ^ ^ E E E E E^f'.i4:;W' (2) 

"=1j1=Ji J2=J2 J3=J3 ^1=0 h=0 l3=0 

where, n(-'^(x) is given by 

2^1 -1 2^2-1 2^3-1 

W = E E E (3) 

Zl=0 l2=0 l3=0 



The index j = (ji, j2, J3) stand for 3-D scales L2/2^^,L2/2^^, and index 1 — (/i, I2, 13) 

for the position of cell < xi < {h + l)Li/2^\ l2L2/2^^ < X2 < {h + l)-^2/2^■^ 

hLs/2^^<Xs < (/3 + l)V2^'- 

In eq.(3), e? 1 is called scaling function coefficient (SFC) of n^(x). They are given by 

ell = / n^(x)(/)j,i(x)o?x = ^ Wi(/)j,i(xi), (4) 

i=i 

where we have used equation (1) in the last step. The scahng functions 0j,i(x) play the 
role of window functions for the cell with volume {Li/2^^) x {L2/2^'^) x (^3/2-'^), and 
located at position 1. Therefore, the SFC e? j is the mean of field n^(x) in the volume 
(Li/2^1) X (L2/2^'2) X {L-i/2^^) at position I. Thus, the term n'^\x) of eqs.(2) and (3) 
actually is the smoothed n'^^x.) by window on scale j. 

The 3-D scaling functions 0j,i(x) can be constructed by a direct product of 1-D scahng 
functions, i.e. 

= (Pjl,hiXl)(pj2,l2i^2)(pj3,l3{X3)- (5) 

In eq.(2), ■^^^■'(x) are wavelets. An 1-D wavelet ipj,i{x) is the modes used to extract 
the fluctuations of a 1-D fleld n{x) on scale L/2^ at position 1. The 3-D wavelets 
are given by mixed direct products of 1-D scaling functions (f>j,i{x) and wavelets '4'j,i{x) as 



follows, 

= 'p3[,hi.^Mj'^Mi^2)^f^Mi^3)^j[,h (6) 

= '^j[,h {^Mi'^M i^2)^j'„l3 i^3W2,j2 

= '4'j[,h i^Mr^h {x2)<f>j'^,h (xsW.js 

Accordingly, ^/'^"/(x) with n = 1,2,3 describe 2-D projected fluctuations on a slice, and 
n — 4,5,6 1-D fluctuations along a line. ■0j/j(x) describes the 3-D fluctuations inside the 
cell 1 = {h, I2, h) on a scale /. 

Hereafter, we consider the projection of density fluctuations onto the space spanned by 
V'jJj (x) only. For simphcity, we ignore the upper index n = 7 in the notations of z/'j^j (x) and 
iyj without confusion. Thus, the wavelet function coefficient (WFC) e?! of n^(x) is given by 

Ng 

ef,! = / ^^WV'j.ilxjdx = Y,Wi%lj^,i{^i), (7) 

i=l 

Therefore, the WFC 1 is the amphtude of the fluctuations of the held n^(x) on scale 
(Li/2^1) X (^2/2^2) X (L3/2J3) at position 1. 

The functions 0j,i(x) and iljyi\x.) form a complete and orthonormal basis in space 
of functions deflned in 3-dimensional coordinate space. The decomposition of eq.(2) is 
information loseless. 

The DWT power spectrum of 6 (x) can be estimated by 
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where (...) is the average over ensemble. The first term on the r.h.s. is the normahzed 
power from n^(x), and the second term corrects for the Poisson noise. The normahzation 
factor n(j, 1) is the mean selection function in the mode (j, 1), i.e., the number density of 
galaxies in the mode (j, 1) when galaxy clustering is absent. As has been shown in paper I 
(Fang & Feng 2000), the factor n(j, 1) can be absorbed into the weight factor by 

n(x) ^ n(xi) 

Thus, equation(8) yields 

^ 2-'l-l 2M-1 2J3-1 2-'l-l 2J2-1 2M-1 ^ ^2^(^^^ 

^ oh+h+js 2^ 2^ 2^ [^j'l] ~ o^I+h+h S ^ ^ / n(:x.) 
^ /i=o ^2=0 i3=o ^ ii=o ^2=0 h=Q "V-^^ 

where the WFC ej^i is defined as 

r n^fx.) 1 



2.2. Selection function in the DWT representation 

By definition eq.(l), selection function fl(x) is an observed galaxy distribution if 
galaxy clustering 5(x) is absent. That is, equation (1) requires to decompose an observed 
distribution n^(x) into two parts: the "background" ri(x), and the fluctuations S{x) upon 
the background. 

This decomposition is not easy if the selection function ri(x) is position-dependent, i.e. 
it will mix with the fluctuations to be detected. With the DWT analysis, one is capable of 
performing this decomposition by a scale-by-scale analysis. 

When we study the fluctuations on a scale j, all x-dependencies of n^(x) on scales 
larger than this scale play the role as a background. Thus, in terms of the scale j, the 
background is given by a smooth of n^(x) on the scale j, i.e. it does not contain information 
of fluctuations on scales equal to or less than j. 
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This background has already be found from equations (2) and (3). Because of the 
orthogonal relation 

j 0j,i(x)^j/,i/(x)cix = 0, if one of j- (i = 1, 2, 3) satisfies j- > ji, (12) 

the function ?t,(-')(x) does not contain any information of fluctuations of modes (j',1) with 
j'i > ji- Thus, to detect the fluctuation power on the scale j, the function n^-'^(x) is 
recognized as a background field. One can identify the selection function as n^-'^(x), i.e. 

2^1-1 2^2-1 2^3-1 

X) = E E E (13) 

Zi=o /2=o ;3=o 



where the superscript j means that this "selection function" is only for the scale j. 

In the plane parallel approximation, selection function depends on ^3 only, i.e. the 
coordinate in redshift direction. From equation (13), the selection function is given by 

n{xs) = n(°°^3)(^3) ^ y: ego,-,^0.3,K^3). (14) 

/3=0 

Prom equation (6), the SFC eoo^y is actually given by an average of n^(x) over the plane 
{xi,X2), and decomposition along X3 direction. 

Using equation (14), the WFC ej^i is now calculated by 

gj,I = E ^(00.3)(^3j ^-^j.'(^^)- ^^^^ 

which presents a simple algorithm for deriving the selection function from observed galaxy 
samples. 

To test this algorithm, we produce mock galaxy samples using N-body simulation, of 
which the details will be given in §4. For simplicity, the selection effect is applied along 
one axis (e.g., direction) of 3-dimcnsional Cartesian coordinates under the plane parallel 
approximation. For the simulation sample in the cubic box with a side length of 256 h~^ 
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Mpc, we replicate the sample along direction, and choose the mock galaxies located 
between 100-356 h~^ Mpc. A selection function is given by 

^(x) = ; l—TT, (16) 

where L is the size of the sample. To be comparable with, for example, the LCRS selection 
function, the parameters are adopted to be L = 500 h~^ Mpc, a = 30 and 6 = 3. 

Fig. 1 displays the DWT diagonal power spectrum Pjjj of particle distributions for 
three typical models with the selection function eq.(14), in which we take js = 7. It shows 
that the power spectrum estimator eqs.(lO) and (15) can perfectly recover the DWT power 
spectrum regardless the selection functions (16). This result will keep unchanged if > 7. 
It means that n^^^'^{'Xs) gives a proper estimation of selection function if is high enough. 
Practically, one can find a properly recovered power spectrum by checking whether it is 
insensitive to j^. 

It should be pointed out that with the developed algorithm the power spectrum eq.(8) 
is normalized scale-by-scale. The fluctuation amplitudes i ( WFCs) on scale j is normalized 
by n(j,l) or n'^^'^'^^^x^) [eq.(15)], which contains fluctuations on all scales larger than j. 
Therefore, the normalization factor generally is scale- dependent. This is different from 
conventional normalization, which is scale-independent. 

If the field is Gaussian, there is no correlation between ejj and n(j,l) or n^^'°'^^\xs)- 
Eq.(8) yields 



It has been shown with the so-called "partition of unity" of wavelets that (l/n^(j,l)) is 
approximately independent of j. In this case, the power spectrum eq.(8) is the same as that 
given by conventional normalization (Jamkhedkar, Bi & Fang 2001). 

However, for non-Gaussian field, especially, if the perturbations on different scales 
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are correlated, the small scale fluctuations given by the WFC e? j generally are correlated 
with large scale fluctuation contained in n(j, 1) or n^^'^'^^\x^). In this case, the power 
spectrum eq.(8) is very different from usual power spectrum which doesn't not sensitive to 
the correlation between perturbations on different scales. 

The power spectrum defined by eq. (8) will benefit to calculate the effect of selection 
function upon the redshift distortion (§3.3). 



3. Redshift distortion in the DWT representation 

3.1. Velocity field 

The redshift distortion is due to pecuhar motion of galaxies. For a given mass field 
5(x), the galaxy velocity v(x) is a random field with mean 

V(x) = (v(x))„ (18) 

where (..)„ denotes ensemble average for velocity fields. V(x) is the bulk velocity at x. 
In linear regime, the bulk velocity is related to the density contrast by 

<^(x) = -T^V-V(x), (19) 

where (3 ~ is the redshift distortion parameter at present, i.e. redshift 2; ~ 0. 

The rms deviation of velocity v(x) from the bulk velocity V(x) is 

{M^)-V,{^)]\ = a%{^). (20) 

In the DWT representation, one can calculate crp„(x) by the variance of the WFCs of 
velocity field, i.e. 1 — J v(x)'0j_i(x)(ix, which is actually the DWT pairwise peculiar 
velocity (paper III). 
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In the paper III, we analyzed the PVD of the velocity fields given by the N-body 
simulation for the CDM family of models. On large scales, the velocity field is basically 
gaussian. It is completely described by its mean V(x) and variance apy{'x). On small scales, 
the velocity field is significantly non-gaussian. The one point function of ej j is lognormal. 
The two-point pairwise velocity correlation functions are non-zero on scales < 1 h~^ Mpc. 
Moreover, on these scales, the pairwise velocities are also correlated significantly with 
density fiuctuations. 



3.2. The DWT power spectrum in redshift space 

The position of galaxy i in redshift space is given by Si = Xi -|- fi>j.(xi)/ifo, where Vr is 
the radial component of v(x). The number density distribution in redshift space is then 



n 



'(s) = Y^^wM^ - - rVr{^)/Ho] = ri^(s)[l + 5^(s)], (21) 



i=l 



where n^{s) is the selection function in redshift space. Similar to equation (8), the power 
spectrum in redshift space is 

^ 2^1-1232-12^3-1 2-'i-l 2^2-1 2^3-1 ^2^(^^) 

" 2Ji+^2+i3 S 5Z 5Z - ^h+h+h ^ ^ ^ / ^fer^^' ^^^^ 

^ /l=0 l2=0 h=0 ll=0 l2=0 l3=0 



where 



''^^'Vj,i(s)rfs. (23) 



The selection function n^{s) can be determined from the observed distribution n^{s) by 
equation (13) or (14). It can then be absorbed into the weight Wi as equation (11) or (15). 
The efi^ect of the diff'crcncc between n^[s) and n(s) will be studied in next section. In this 
section, we will ignored this effect. Thus, using an auxiliary vector J, equations (21) and 
(23) yield 



J 



f 1 = Y:^ifds6n{s-^, + tVj) e^-^--'^('^*)/^°^j,i(s) |^^^ , (24) 

i=l 



where Vj is gradient operator on J. 



Subjecting equation (24) to an average over ensemble of velocity, if the velocity field is 
gaussian, we have 



j=o 



i=i 
i=i 

For simplicity, we will use e?i for (e?i)„ below without causing confusion. 
If we consider only linear effect of the bulk velocity, equation (24) gives 



(25) 



j=o 



1 + i—i ■ rV;(xi) 

-no 



ef 1 = ^ / (is (s - Xj + iVj) 
1=1 ■' 

= Y.^^ dsV^j,i(s)e-(^/2)'^-(^)(^-^)'5,,(s - X,) 

i=l ■' 
1 f 

- TrlZ'^i dsr-[VsVr{s + iVj)SD{s-Xi + iVj)]e~ 
^0 i=i J 



-(l/2)a2„(xi)(J-f)^ 



V'j.Ks) 



J=0 



(26) 



(l/2)a^„(xi)(J-f)%, 



3=0 



Neglecting the terms of the order of K(x)5(x), and using the linear relation eq.(19), 
equation (26) gives 



+ 



J dsV^j,i(s)e(V2)^|.(s)(r.V)2^g(g) (27) 

p J (isV^j,i(s)(r- V,)V-2e-(i/2)a^„(s)(f.v)2^s(s) 



Because all operators in the integrand of equation (27) are nearly diagonal in the DWT 
representation (Farge et al 1996), equation (27) can be rewritten as 



where 



el = (1 + (^S,)s?%, 



Si- /V'j,i(x)(f-V)V-Vj,i(x)ci: 



X. 



(28) 



(29) 
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and 

The method of calculating Sj and sj^ is presented in Appendix A. 

Substituting equation (28) into equation (22), we have the redshift distorted power 
spectrum as 

pf={i+psirsrpy (31) 

where Sf^ — [s?^]^. Above equation formulates the redshift distortion effect in DWT 
expression. This is the basic formula for our redshift distortion analysis. Usually, the factor 
(1 + /5Sj)^ is called linear redshift distortion, and S^"^ called non-linear redshift distortion 
due to the pairwise velocity dispersion. However in our derivation, the two parts are not 
treated separately. 

This derivation can be generalized to any velocity fields, which are not simply described 
by equations (18) and (20). For instance, if the pairwise velocities are correlated, i.e. 

(b,(x) - y,(x)]h(x') - T/,(x')]). = <7j,(x - x'), (32) 

the pairwise velocity dispersion factor becomes 

Sf^^JJd^d^' (33) 

^/^j,l(x)^/^j,l(xOe[(V2)-?Jx)(^V)2+a|Jx-xO{^V)(^VO+(l/2)a?JxO(^VO^]^.j(x)^/;j^ 

where V is gradient operator on x'. 

3.3. Effect of selection functions 

The theory of redshift distortion presented in last section did not consider the effects 
of selection functions, n^{s) and n(s). Since selection function relies on the radial distance. 
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it could be also a source of the anisotropy of power spectrum with respect to redshift 
direction. It should be taken into account when analyze observed samples. 

In the linear approximation of Vr, equation (21) gives 



n 



^« 1, _ 1 



s) = n^(s) - J2 Tr^-(^i)^ • - X,) = n^(s) - — r • VK(s)K(s)]. (34) 



=1 -no -no 



In the second steps, Vr is replaced by Vr, as only the bulk velocity is considered. Using 
equations (1) and (21), equation (34) yields 

where the second order term (5(s)l/r(s) is ignored. 

The term -1 on the r.h.s. of equation (35) does not contribute to power spectrum, 
because of / tlj^^i{x)dx — 0. In the hnear approximation, n(s)/n^{s) ~ 1 + 0{vr) and 
therefore, the factor n^(s)/n'^(s) in the third term of the r.h.s. of equation (34) can be 
approximated as 1. 

The second term in the r.h.s. of equation (35) contains a linear term of Vj., i.e. the 
same order as the third term. However, with the selection function equation (13), we have 

j {x)ipj^i{x)dx = j n^^ \x)ipj^i{x)dx = if one of ji {i = 1, 2, 3) satisfies ji > j[. (36) 

For plane parallel approximation, j' = (0, 0,j3), the above equation is always hold if we 
study fluctuation powers on scales ji > or j2 > 0. Actually, in this case, n^(x) or n*^j'^(x) 
depend only on xs, and thus their projections onto bases ipj^i{xi) or ipj,i{x2) {j > 0) are 
always null. Similarly the selection function in redshift space is also only dependent on x^. 
Accordingly, in the plane parallel approximation, we have 

/fi(s) 
-^V^j,ids = ifiiori2>0. (37) 

which implies that the second term in the r.h.s. of equation (34) also has not contribution 
to the projection on bases ■0j,i(^i) or '0j,i(^2)- 
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In fact, the condition of plane parallel approximation is not necessary. For radial 
redshift, the selection function is still a 1-D function. Its projection onto the bases V'j,! on 
the celestial spherical surface is still zero. Thus, the linear redshift distortion mapping is 
given by 



S^{s) ~ 5{s) - 



1 



Hon{s) 



V[n{s)vr{s)]. 



(38) 



Hereafter we will use x for the variable s. It will not cause confusion as the superscript S 
stands for redshift space. 



In the plane parallel approximation, equation (38) becomes 



5^(x) = 5(x) 



1 dvs 1 dlnn{x3) 



Ho 8x2, Ho dx3 
Using equation (19), vs can be represented by 5(x), we have then 



(39) 



^ ^ dlnn{x3) d 
dx'i dx:i dx:t 



5(x). 



(40) 



The differential operator of the second term in the bracket of equation (40) is nearly 
diagonal in the DWT representation (Farge et al 1996.) Thus, we have 

dlnn^Xs) 



e* = ,^(l + /55j)ej,, + /3- 



dxs 



(41) 



where dlnn (x^) / dx^l^^i means the value of dlnn {x3)/dxs in the mode The coefficient 



Qj,i3-i'^ are defined by 

/O 

The calculations of Qj,is-i' are given in Appendix A. 



(42) 



Since Qj^o — (Appendix A), the first and second terms on the r.h.s. of equation (40) 
are not correlated, we have 

n 2 

d\nn{x3) 



dxs 



(43) 
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For uniform fields, (|e'^j_iP) and (|ej_ip) are 1-independent. Hence, equation (43) gives the 

relation between the DWT power spectra in redshift Pf and real spaces Pj as 

_ n 2 

din 77,(0:3) 



[1 + PSiY + 



dxs 



j.i. 



> Pi- 



rn 



which quantifies both redshift distortion and selection effect on the DWT power spectrum. 



Using inequality equation (A15), we can show that if 

din 71(2:3) 

^3 



we have 



P 



dx3 



j.i. 



(45) 



(46) 



which states that the selection function term in equation (44) is even less than the second 
order terms P^S? if the selection function varies with X3 slowly. In this case, the selection 
function does not significantly disturb the linear redshift distortion described by equation 
(31) when we perform a DWT power spectrum analysis using equations (10) and (15). 

The Poisson noise term in equation (10) is not affected by the redshift distortion, i.e. 
it is the same as the Poisson noise term in equation (22). Because both Poisson noise terms 

of equations (10) and (22) linearly depend on 7z(x) or 77'^(x), for an ensemble average, we 
have (n(x)) = {n^{x)). 



4. Simulation samples 

To demonstrate the redshift distortion in the DWT representation, we use the N-body 
simulation samples like that in Paper II. The model parameters used are listed in table. 1. 
We use modified AP'^M code (Couchman, 1991) to evolve 128^ cold dark matter particles 
in a periodic cube of side length L. The linear power spectrum is using the fitting formula 
given in Bardeen et al. (1986). 
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In our simulation, we use the so-called "glass" configuration to generate the unperturbed 
uniform distribution of particles, and the Zel'dovich approximation to set up the initial 
perturbation. The triangular-shaped cloud (TSC) method is used for the mass assignment 
on the grid and the calculation of the force on a given particle from interpolation of the 
grid values. We take 600 total integration steps from zi — 15 for the SCDM model, and 
Zi = 25 for ACDM and rCDM down to 2; = 0. The force softening parameter rj in the 
comoving system decreases with time as 77 oc l/a{t). Its initial value is taken to be 15%, 
and the minimum value to be 5% of the grid size, respectively. 

5. /3-estimators and reconstruction 

5.1. A test of the linear approximation 

We try to estimate the redshift distortion parameter /3 using equation (30). Obviously, 
the precision of the (5 determination is dependent on the linear relation of equation (18), 
which is valid only on the scales where the bulk velocity can be described by linear or 
quasi-linear density perturbations. It has been already realized that the non-linearity of 
the relation between (5(x) and V(x) will be significant on small scales (e.g. Kudlicki et al. 

Table 1 

Model L/h~^Mpc ^ T as /3 realizations 

SCDM 256 1.0 0.0 0.5 0.55 1.0 10 

rCDM 256 1.0 0.0 0.25 0.55 1.0 10 

ACDMl 256 0.3 0.7 0.21 0.85 0.49 10 

ACDM2 480 0.3 0.7 0.21 0.95 0.49 6 
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2000). Therefore, it is necessary to have an estimation of uncertainties due to the non-hnear 
5(x) — V(x) relation. 

The non-hnear effect can be estimated by a "blueshift" (change the sign of velocity 
from v(x) to — v(x)) distorted power spectrum. Prom the derivation of equation (31) in 
last section, it is easy to show that, the blueshifted power spectrum is 

Pf^il-PS^rsf^^Py (47) 

Which has the same non-linear redshift distortion effect Sf^ as the redshift power spectrum. 
The difference between the redshift [equation (31)] and blueshift [equation (47)] distorted 
power spectra is a sign of the linear term /3. Combining equations (31) and (47), one can 
determine (3 as 

' - ^ WW * ^ 

We test equation (48) by simulation samples in the ACDM model. First, we calculate 
the diagonal power spectra of particles, P^^j and P^^p in the plane parallel approximation. 
We then find /3 by equation (48). The result is presented in Fig. 2. The density parameter 
used for the simulation is Q = 0.3, or P — 0.49. Fig. 2 shows that the values of /3 given 
by equation (48) is generally overestimated. But the overestimations are no more than 
10% on scales k < 0.5 h Mpc~^, and about 20% for k > 1 h Mpc^^. It should be pointed 
out that the goodness of the resulted /? on small scales k > 1 h Mpc~^ doesn't mean that 
the linear relation eq.(19) can be used to describe the redshift distortion on these scales. 
The non-hnear effect on these scales can not be ignored, and the redshift distortion is 
dominated by the term <Sj^^. The goodness for A; > 1 h Mpc~^ shown in Fig. 2 is due to 
that the non-hnear effect is significantly repressed by the test eq.(48) based on an assumed 
blueshifts. 
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5.2. Symmetry of the DWT quantities 

At the first glance, it seems impossible to estimate (3 using the observed redshift 
distorted power spectrum P-^ only, as all the quantities such as Pj, and S^^ (or a^^) 
included on the r.h.s. of equation (30) are unknown. Usually, equation (30) is used to fit 
the observed power spectrum with model-predicted i^, and then determine the parameter 
/3 from the most-likely-fitting. 

But we try to search for /3-estimators which depend on model as less as possible. To 
achieve this, we take the advantage of the wavelet analysis: the DWT modes of V'jij-aj'sC^) 
are not rotational invariant, but cyclic permutational invariant. As a consequence, all the 
DWT quantities in the redshift distortion equation (30), i^, P-^ , 5j, and "^j^^, are dependent 
on the three indexes (ji, j2, Js), rather than the length scale of mode j only. The quantities 
Pj, ij"^, 5j and S^^ satisfy the following symmetry. 

1. If cosmic density and velocity fields are statistically isotropic, the DWT power 
spectrum in real space is invariant with respect to the cyclic permutations of index 

j = (ii,i2,i3), i-e. 

2. In the plane-parallel approximation, e.g., coordinate is in the redshift direction, 
we have 

pS _ pS /irri\ 

jl,h,j3 j2jlj3 / 

^31,32,33 ~ ^32,31,33 i^^) 

qPV _ qPV (r<2) 
h,32,33 32,31,33 \^ ) 

3. If (Tp^ is constant, i.e. scale- independent, following equation (AlO), we have 
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i.e. Sf^ is independent of ji and j2. More generally, equation (53) also hold even when the 
radial correlation of pairwise peculiar velocity is considered [equation (33)]. 

Equations (49) - (53) provide the base of designing the /3-estimators with the DWT 
power spectrum. 



5.3. /3-estimator with scale- independent cr. 



pv 



Assuming that the pairwise velocity dispersion is scale- independent, equations (31) 
and (49) give 



^33233 ^ (1 + f^Sjj^j^) 

^]%32 i'^+^^3j2.hy 



-I 2 



qPV 

33233 
CPV 
33332 . 



(54) 



For a given pair (j2, Js), ^fj^jz/^fhh ^ constant, and thus the r.h.s. of equation (54) 
depends only on two parameters (3 and Sfy.lSfy.. Accordingly, these parameters can be 
found by the best fitting of the r.h.s. of equation (54) with observed ratios -P;§2i3/-^jj3i2' 
J = 2,3... 

A numerical example of this fitting is demonstrated in Fig. 3, in which we take 
(j2,j3) = (2,3), and j = 2. ..7. We calculate Pjj^js/ Pjj^ji ACDM simulation samples, 

and the best fitting yields the values of /3 = 0.53 ± 0.25 and Sf^^/Sf^ = 0.80 ± 0.04. The 
precision of the estimator eq.(53) probably is not better than 20%. This shows that the 
assumption of a constant cxp^ is not too bad, but its effect on the /3-determination cannot 
be neglected. 



5.4. /3-Estimator with scale-dependent Gpy 

If (7py is scale-dependent, equation (54) is no longer correct in general. In this case 
S'j^j^ja depends on the three index (ji, j2, Ja), rather than j'3 only. So the relation eq.(30) 
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is obviously not enough to extract (5 from measuring redshift distorted power spectrum P.^ 
only. To search for an appropriate algorithm for (5 estimation, we shall first consider the 
property of S^^ . 

At first, the pairwise velocity dispersion factor S^^ is determined by an isotropic 
function (jp^(x), which is scale-dependent. Let us consider the modes of LjV^ < 
where j± is defined by 2^-?-^ = 2^-^i -|-2^-'2. In this case, the scales of these modes are dominant 
by ja, and therefore, (Jp„(x) is also dominant by js- That is, for a given js, Sf^ will keep 
constant if j± > j^. 

To test this expectation, we calculate Sf^ in the ACDM model. The result is shown in 
Fig. 4. One can see from Fig. 4 that S^^ is generally dependent on indexes j± as well as j^. 
However, for a given js, Sf^ almost keeps constant in the range of j± > j^. Since j = 2, 3 is 
the two largest scales of the samples, 'S'j^]^2 ^hJ23 ^^^P constant very well for all ji, j2. 



Using this result, we can design a /9-estimator as follows 

Pf'23Pf32 {l + PSy23m + PSjs,r 



or 

(3 ~ 



' pS pS \ 1/2 
^j23^j'32 \ _ ^ 

Pr23Pj32) 



(55) 



(56) 



Obviously, the first factor on the r.h.s. of equation (56) is to measure the difference between 
spectra with the same scales, but different shapes of the modes. The estimator eq.(56) is 
model-free, as equation (56) is based only on property that two DWT modes of (j±, J3) and 
(ii) is) have the same scale, but different shapes if > is and j'^ > j^. 

To apply the estimator (56), we take i = 2 and / = 7, because modes with i = 2 and 
j' = 7 have largest difference in the shape, but very small difference in the scale. This 
estimator yields the values of /? = 0.47 ± 0.18 for ACDM sample with simulation box L — 
480 h-i Mpc; 0.93 ± 0.22 for SCDM sample with L = 256 h-^ Mpc; and 1.00 ± 0.34 for 
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rCDM sample with L — 256 h ^ Mpc. Considering the uncertainty of non-hnearity is about 
20% (§5.1), the estimator eq. (56) works very well. 



5.5. Estimation of 

After the j3 estimation, we can also estimate cip^ on scale J3 using off-diagonal power 
spectrum in redshift space. From §5.4, we know that S^^ is almost constant in the range 
of 3± > h- So we have Sf^^ = Sfy and Sf^ = S^^. 

We can calculate the ratio of pairwise velocity dispersion factor Sf^^ /Sf^ by 

Pf^jjl + pSjrrY ^ ^ Sj^^ ^^^^ 

On the largest scale J = 2, the nonlinear redshift distortion due to the pairwise velocity 

dispersion is negligible, we have 5'^2 — 1- Thus, we have all the diagonal members of the 

pairwise velocity dispersion factor -S'^J. The parameter on the scale j can then be found 

by = with equation (AlO), i.e. 

-1 00 

C.3 = ^ E mns/2^^)\'eM-im4vi'^^ns/L,n (58) 

n3=oo 

For the ACDM simulation sample, S^'^ is shown in Fig. 5. The on scales j = 3. ..7 
is shown in Fig. 6. Although the values of apy shown in Fig. 6 are correct in average, but it 
does not match with the direct measurement of apy given in Paper III. Especially, on small 
scales j — 6 and 7, the values of ap^ are significantly lower than the direct measurement. 
This is not unexpected. The factor equation (58) is obtained under the assumption of 
gaussian velocity field (§3.1). However, we have shown in Paper III that the velocity field is 
actually intermittent on small scales. The PDF of ej j is not gaussian, but lognormal. 

Nevertheless, the factor ^j^^ is still good for the redshift-real mapping and /9-estimation, 
because in these calculations, we used only the symmetric properties and scale-dependence 
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of Sf'^, but not the details of ap^. This point can also be seen from reconstruction of power 
spectrum in real space from that in redshift space. Using the P estimated by equation (56), 
and Sj'^^ by equation (57), one can reconstruct the diagonal DWT power spectrum in real 
space Pjjj from P^jj through equation (31). Fig. (7) compares the recovery of Pjjj with the 
original diagonal DWT power spectrum in the ACDM model. The result shows that the 
algorithm of reconstruction is reliable. 



6. Conclusion 

We established the mapping between the DWT power spectra in real and redshift 
spaces. Prom equations (31) and (44), the mapping in the plane parallel approximation is 

n 2 

d\nn{x3) 



E Qh.-. S^Py (59) 



dX2 

which includes the effects of (1) bulk velocity (the term of Sj), and (2) selection function 
(the term of n{x3)) as well as (3) pairwise peculiar velocity (the term of Sf^). 

In the Fourier representation, the redshift distortion mapping is axially symmetric with 
respect to the redshift direction. This is because the Fourier mode is rotational invariant, 
and the redshift distortion produces an anisotropy in the spatial distribution of galaxies 
between the line of sight and directions other than it. An isotropic statistics plus an axially 
symmetric violation will result in an axial symmetry. On the other hand, all redshift 
distortion factors in equation (59) are no longer axial symmetric, which is due to the 
modes j being not invariant under rotational transformation, but for cyclical permutation. 
Therefore, the symmetry of the DWT redshift distortion results from a cyclical permutation 
plus an axially symmetric violation. 

By virtue of this feature, we develop /? estimators, which are mainly based on the shape- 
dependence of the redshift distorted DWT power spectrum, but not the scale-dependence. 
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The estimators eqs.(54) and (56) look similar to the quadrupole-to-monopole ratio method 
based on the Fourier power spectrum in redshift space, which gives a model- free estimation 
of /? (e.g. Cole et al, 1994, Fisher et al, 1994, Hatton & Cole, 1998, Szalay et al, 1998). 
However, the estimators eqs.(54) and (56) is different from the quadrupole-to-monopole 
ratio method. The latter considered only the hnear redshift distortion effect, but not the 
non- linear redshift distortion effect. If considering the non-linear redshift distortion effect, 
the quadrupole-to-monopole ratio method needs a fitting of observed redshift distorted 
power spectrum with model-predicted power spectrum. On the other hand, the estimators 
eqs.(54) and (56) considered both linear (/3Sj) and non-linear (Sf^) redshift distortion 
effects. These /3-estimators are model-free even when the non-linear redshift distortion 
effect is not negligible. We test the /3-estimators using N-body simulation samples. The 
result shows that regardless the pairwise velocity dispersion is scale-dependent or not, 
the /9-estimators can yield the correct number of /? with about 20% uncertainty. We also 
develop an algorithm for reconstruction of the power spectrum in real space from the 
redshift distorted power spectrum. The numerical tests also show that the real power 
spectrum can be reasonably recovered from the redshift distorted power spectrum. 

LLF and YQC acknowledges support from the National Science Foundation of China 
(NSFC) and National Key Basic Research Science Foundation. 



A. Calculations of Sj, s^"", and Qj.i-i' 



Let us consider plane-parallel approximation, i.e. coordinate is in the redshift 
(2;)-direction. The linear redshift distortion Sj equation (28) gives 

^2 



Because 1-D wavelet ipj,i{x) is given by dilating and translating the basic wavelet ip{r]) as 

(A2) 

The Fourier transform of tpj^i is 

= 7 E (A3) 



4vW = (4V^'v^W2^>"^'™''"'' (A4) 



and 

lit ■ 1 in 1 = 1 

.23 J 

where ip{n) is the Fourier transform of the basic wavelet 

ijj(n) = / ^/;(77)e-^''^"''d77. (A5) 

— oo 

The function |'^(n)p is shown in Fig. 1 of Yang et al. (2001). 
Thus, equation (Al) becomes 

-^^1^2,^3 = '^j^Th+h ^^^^ 



\2 



ni,n2,n3=oo + ("-2/^2)^ + [n^l L^Y 

Since ■^(n) is non-zero only around two peaks at n = ±np, the summation of equation (A6) 
actually only over few number of ni around ±np2^'. 

If Li — L2 — — L, equation (A6) becomes 

1 °° 772 . 

^.....3 = E ,2 , J , ,2 |^(ni/2^-^)^(n,/2^-)^(n3/2^3) |2_ ^^7) 
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For diagonal modes, i.e. ji — j2 = J3 = j we have 



Sj,j,j = I (A8) 



A.2. sf 



If (T^(x) is independent of x, s?^ [equation (28)] can be calculated by 



-1 00 

= ^^IT^ E |ViK/2^0Vi(ri2/2^-^)^(n3/2^-^)rexp[-(l/2)a^,(r • nf], (A9) 

ni,n2,n3=oo 

where vector n = 27r(ni/Li, ^2/^2, ris/Ls). In plane-parallel approximation, we have 

1 00 

C.3 = ^ E mns/2^^)\'eM-im4vi'^^ns/L,n (AlO) 

n3=cx) 

The summation of equations. (A9) and (AlO) also runs only over few number of rii around 



A. 3. Q\,i-v 



Similarly, for given by [equation (41)] 



/O 
i'},h,l2,l3 i^)-Q^'^~^^i,h,l2,l's W^X' (All) 



we have 



ni,n^3=oo 27r[(ni/Li)^ + (na/La)^ + (na/La)^] 
Equation. (A12) gives 

gj,i3_. = 0, if - ^3 = 0. (A13) 
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Prom equation (A12), we have 

Q\ 1-1' = . }. ^. (A14) 

ni,n^;r3=oo ^3 27r[(ni/Li)2 + (na/La)^ + (713/^3)2] 

~ 27rnp2-?3 2-?i+J2+-?3 

~ (n3/L3)2sin[27rn3(/3 - /')/2-'-^l , , . x - . , . ^,9 

for the last step, we consider that ip{n3/2^^) requires ~ 2^^np, and rZp ~ 1 being the peak 
of ^(n). Since J2f=o^ sm{27ml/2^) < 1, equations. (A7) and (A14) yield 
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Fig. 1. — The diagonal DWT power spectrum (square) measured in simulation samples for 
models SCDM (upper panel), rCDM (central panel) and ACDMl (lower panel). The error 
bars are 1-cr variance from 10 realizations for each model. The solid lines are the original 
power spectrum 

Fig. 2. — The values of (3 estimated by the redshift and blueshift power spectra in the 
ACDM2 simulation samples. The error bars are 1-a variance from 6 reahzations. 

Fig. 3. — Pjmzl ■^jjzj2 3 ACDM2 simulation samples. The error bars are 1-cr variance 
from 6 reahzations. The sohd hne is given by a best fitting with eq.(52) 

Fig. 4. — S^^ vs. j± of ACDM2 simulation samples. The error bars are 1-a variance from 
6 realizations. 

Fig. 5. — S^jj /S2-^ of ACDM2 simulation samples. The error bars are 1-a variance from 6 
realizations. 

Fig. 6. — PVD estimated from the ACDM2 simulation samples. The error bars are 1-a 
variance from 6 realizations for each model. 

Fig. 7. — Reconstructed DWT power spectrum for the ACDM2 simulation samples. The 
error bars are 1-a variance from 6 realizations. 
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